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We extend the method of antisymmetrized molecular dynamics to investigate dineutron corre¬ 
lation. We apply this method to ^°Be as an example and investigate the motion of two neutrons 
around a largely deformed ®Be core by analyzing the two-neutron overlap function around the core. 

We show that the core structure plays an important role in dineutron formation and expansion from 
the core and that the present framework is effective for the studies of dineutron correlation. 


I. INTRODUCTION 


Many exotic phenomena have been found in the 
neutron-rich nuclei, and more have been suggested by 
both theoretical and experimental studies. One such phe¬ 
nomenon is dineutron correlation, a strong spatial cor¬ 
relation between two neutrons coupled to a spin singlet. 
Although the two neutrons are not bound in a free space, 
strong dineutron correlation has been theoretically sug- 


S sted, e.g., in a low-density region of nuclear matter 
,12 or in the neutron-halo or -skin regions of neutron- 
rich nuclei [IS3- These studies clarify that the strength 
of dineutron correlation reflected in the dineutron size 
changes significantly depending on circumstances such 
as nuclear density and potential from the core. In addi¬ 
tion, the dineutron and diproton correlations have been 
intensively discussed via 2n and 2p emissions from the 
unbound nuclei in connection with recent experiments 


Preceding studies have investigated dineutron correla¬ 
tion in the ground and excited states in certain nuclei, but 
they have not conducted a systematic investigation, and 
the formation mechanism underlying dineutron correla¬ 
tion and the dynamics of one or more dineutrons are not 
well understood. To investigate in detail the dineutron 
motion with resMct to the core, a core-|-2n three-body 
model is useful Miiini- However, in three-body 
models, it is somewhat difficult to take various structure 
changes into account (including excitation, deformation, 
and clustering), and inert, spherical cores have been as¬ 
sumed in most of the preceding studies. Core excitation 
and deformation can affect dineutron correlation, and it 
is necessary to consider changes in the core structure ex¬ 
plicitly for a systematic investigation of dineutron corre¬ 
lation. 


In our previous studies, using dineutron condensate 
(DC) wave functions [TsIhtII . we showed that dineu¬ 
trons in nuclei are fragile and easily broken. To study 
the properties of dineutron correlation, the following 
dineutron-breaking effects should be considered. First, 
the dissociation of a spin-singlet two-neutron pair due 
to the spin-orbit potential from the core, as discussed 


in Refs. [T3, HO, HH. Without the spin-orbit potential, 
two neutrons tend to be coupled to a spin-singlet pair 
because of the spin-singlet s-wave attraction. However, 
in reality, two valence neutrons at the nuclear surface 
feel the spin-orbit potential from the core and tend to 
occupy the LS”-favored orbits, resulting in mixing of the 
spin-triplet pair, i.e., reduction of the spin-singlet dineu¬ 
tron component. We call this effect at the surface due to 
the spin-orbit potential the “AS' dissociation” of a dineu¬ 
tron. Second, the swell in the size of a dineutron in the 
region far from the core. Two neutrons are not bound in 
a free space; thus, dineutron correlation vanishes in the 
asymptotic region far from the core, meaning that the 
dineutron size becomes infinitely large as the distance 
from the core increases, as discussed in Ref. [HI- We 
call such a breaking effect at the region far from the core 
“dispersion” of a dineutron. 

To clarify the properties of the dineutron correlation in 
neutron-rich nuclei, we use the extended method of an¬ 
tisymmetrized molecular dynamics (AMD) The 

AMD framework can describe various structures such as 
deformation and clustering in general nuclei, and it is 
suitable to describe various core structures. We extend 
the AMD framework to investigate the degree of dineu¬ 
tron formation at the surface and the degree of expansion 
of a dineutron tail at the farther region while taking into 
account the above-mentioned breaking effects of a dineu¬ 
tron based on the concept of a core-|-2n. In addition, 
to analyze the detailed two-neutron motion, we propose 
a method that enables us to visualize the two-neutron 
spatial distribution around the core. As the first step, 
we apply the extended AMD framework and the method 
of analysis to ^°Be, which has a well-deformed ®Be (2a) 
core, and discuss the effect of core structure change on 
the dineutron correlation in this study. We show that 
the core structure significantly affects the dineutron for¬ 
mation and distribution around the core. 

This paper is organized as follows. In Sec. m we ex¬ 
plain the framework used to describe and analyze the 
dineutron correlation in neutron-rich nuclei. In Sec. m 
we apply the framework to the ^°Be system and discuss 
the dineutron formation and distribution around the ®Be 
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core. We summarize our work in Sec. m 


B. Core+2n AMD wave function 


II. FRAMEWORK 

We consider a total A-nucleon system composed of an 
{A — 2)-nucleon core and two valence neutrons and in¬ 
vestigate the two-neutron motion around the core. We 
describe our framework in this section. 


A. AMD wave function 


First, we explain the AMD wave function. An AMD 
wave function for an A-nucleon system is given by the 
Slater determinant of A single-particle wave functions; 


^AMs{Z) 
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A{ipi ■ ■ ■ ipA} ■ 


( 1 ) 


Here A is an antisymmetrizer and (pi {i = 1,..., A) 
are the single-particle wave functions composed of the 
Gaussian-type spatial part (j)i, the spin part and the 
isospin part n as follows: 


>Pi =4'iXin, ( 2 ) 



Xi = + (4) 

Ti = p or n. (5) 


Z = {Yi,...,y'A,4i,...,lA} in Eq. Q are the varia¬ 
tional parameters that characterize the Gaussian centers 
and spin orientations of the A nucleons, v is the Gaus¬ 
sian width characterizing the size of the single-nucleon 
motion, which is generally common for all nucleons in 
the standard AMD framework. 

In the present AMD framework, the widths can differ 
from one another (z/^ also has the label of the ith single¬ 
particle state) and we choose i^i = Vc {i < A — 2; core) 
and Vi = V 2 n (* > A — l;two valence neutrons), as in 
Ref. (ssj. Note that in the case where all the Gaussian 
widths differ, the center of mass motion cannot be sep¬ 
arated exactly. Therefore, we treat the center of mass 
motion approximately by reducing the center of mass ki¬ 
netic energy from the total energy, as described in IIII Al 


We comment upon the relation between the present 
AMD framework and the similar fermionic molecular dy¬ 
namics (FMD) framework [s^. Our AMD framework, in 
which all the Gaussian widths are not identical, is indeed 
quite similar to the FMD framework; however, we use 
phenomenological effective interactions in the Hamilto¬ 
nian (Eq. (12011 '). whereas, in a recent FMD calculation, 
they used effective interactions derived from realistic in¬ 
teractions via the unitary correlation operator method 


We aim to efficiently describe dineutron correlation in 
neutron-rich nuclei. To this end, we prepare the total 
wave function in two steps. The A-nucleon system is 
regarded as a core composed of A — 2 nucleons and two 
valence neutrons. The two-step treatment of this system 
is explained below. 

First, we prepare the core wave function composed of 
A — 2 nucleons. 

<^Le = A{pt---,V^A-2}- ( 6 ) 

k is the label characterizing the core structure, and 
the core wave function is specified by the parameter 
sets {Y^,..., K^_ 2 ,li,...,^ a- 2 }- We superpose the 
core wave functions consider the core structure 

change. Here, we consider the ®Be core, as shown later. 
We use the 2 q; cluster wave function for the ®Be core for 
simplicity and choose the a-a distance da to parameter¬ 
ize the core structure. We can simply generalize the core 
wave function to the AMD wave function. For example, 
in principle, it is possible to prepare the deformed cores 
by performing the /3-constraint variation, and the defor¬ 
mation parameter /3 is chosen as the parameter for the 
core structure. 

After preparing the basis wave functions of the core, 
we express the A-nucleon total wave functions by dis¬ 
tributing two neutrons around each core wave function: 

4>total = -^A ^A-I'Pa} , (7) 

with being an AMD wave 

function having the parameters 

{Y\, ..., Ya-1, Ya, ..., U-iAa}, 

where the parameters for the core are fixed. We 
perform energy variation only on the parameters 
a-1, Ya,$a-i,^a} for the two valence neutrons 
to minimize the total energy under the condition 
Y, = (0,0,0) with Y, = Y>^ - (Ya-1 + Ya)/2 for 
i < A — 2. The energy variation for two valence neutrons 
is performed under the constraint on the distance of the 
center of mass, (Ya-i + Ya)/2, of two neutrons from 
the core, which we call d 2 n- 


d2n 


i(YA-i + YA) 


1 

A- 2 


E 

i<A-2 


( 8 ) 


Performing the d 2 ri-constraint energy variation for d 2 n = 
d 2 n,m ijn = 1,...), we calculate the two-neutron wave 
functions near and far from the core. The variations of 
Ya -1 and Ya are performed independently so the two 
neutrons near the core tend to have the opposite mo¬ 
mentum (i.e. the imaginary parts are opposite to one 
another) and tend to be broken by the LS dissociation 
effect due to the spin-orbit potential from the core. On 
the other hand, two neutrons far from the core tend to 
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have almost the same position, indicating the expansion 
of a spin-singlet dineutron far from the core owing to the 
energy variation. A similar constraint is proposed for a 
cluster in Ref. [s^ and can efficiently describe a cluster 
development. 

For the valence neutrons, we choose different values 
of the Gaussian width from those of the core nucleons 
{v 2 n ^ Vc). Using the different values of the Gaussian 
widths V 2 m the dineutron dispersion effect is taken into 
account. Moreover, the description of the spatially ex¬ 
panded tail parts of single-particle wave functions of va¬ 
lence neutrons such as the neutron-halo tail should be 
improved using multi-range Gaussians. 

We describe the total system by superposing all wave 
functions with core deformation (fc), core-2n distance 
{d 2 n) and 2n size {v 2 n)'- 

^A{Jtt) = kPM totaxi^ ( 9 ) 

K k d2n ^2n 

Pmk parity and angular momentum projection 

operator to the eigenstate of Jtt. The coefhcients ckk, 
(k is the abbreviation of {fc, d 2 m i^ 2 n\) are determined by 
diagonalizing the Hamiltonian (Eq. (l20l) h Superposing 
the wave functions with the set of {k,d 2 n,V 2 n}, we can 
take the core structure change into account as well as 
the two-neutron motion near and far from the core while 
varying their size. 

C. ^°Be AMD wave function 

Here, we specify the form of the AMD wave function 
for ^°Be. The explicit form of the ^°Be wave function is 
given by 

$iOBe(0+) (10) 

K da d2n ^2n 

^°Be is regarded as a ®Be core with two valence neutrons. 

The ®Be core is described with the 2 q: cluster wave 
function for simplicity. An a cluster is composed of four 
spin- and isospin-saturated nucleons in (Os)"' with the 
same Gaussian center and width. We choose a-a dis¬ 
tance da as the parameter characterizing the core struc¬ 
ture (fc = da in Eq. Q). We prepare three types of cores 
with da = “2, 3, 4 fm to describe the core structure fluc¬ 
tuation. The Gaussian widths of the core nucleons are 
fixed to Vc = 0.235 fm~^, giving the lowest energy of an 
a. 

Then, we distribute two neutrons around each core 
with da = 2, 3, 4 fm under the d 2 n-constraint (Eq. (O). 
Here, we choose d 2 n to be d 2 n = (m = 

l,---,TOmax)- We define = dmax- We set 

(TOmax, rfmin, rfmax) = (8,1,5). We verified that the '°Be 
ground state energy is lowered by only 100 keV at most 
if the number of bases or the maximum value of dmax is 
increased. When d 2 n is small and two neutrons are dis¬ 
tributed near the core, the LS dissociation for a dineu¬ 
tron plays an important role and the two neutrons mainly 


occupy the lowest allowed orbit 0 p 3/2 to gain the spin- 
orbit potential from the core. On the other hand, when 
d 2 n is large, the tail of the dineutron center of mass mo¬ 
tion far from the core can be well described. 

In each d 2 n-constraint variation, we fix the Gaus¬ 
sian widths of the two valence neutrons to be V 2 n = 
0.235, 0.125, or 0.08 fm“^ and perform variation for each 
V 2 n value. The superposition of three types of bases with 
different U 2 n values for each d 2 n can describe the dineu¬ 
tron size change. 

At the nuclear surface, the LS dissociation of a dineu¬ 
tron is predominant and dineutron correlation is sup¬ 
pressed. To more effectively describe the dineutron for¬ 
mation at the surface, we additionally superpose the 
bases prepared with the d 2 n-constraint variation without 
the spin-orbit force {vls = 0 MeV) when d 2 n is small. 
Here, we prepare the bases without the spin-orbit force 
when c? 2 n < 2 fm (four bases) for each {da^ V 2 n) set. 

We summarize the number of bases used to describe 
'°Be. We use three types of ®Be cores {da = 2, 3, 4 fm) 
and three widths of the two valence neutrons for each core 
(i^ 2 n = 0.235, 0.125, 0.08fm“^). We perform the energy 
variation on the two valence neutrons in each {da,V 2 n) 
set under the d 2 n constraint. We choose eight values for 
d 2 n with Vls = 1600 MeV (this value will be explained 
later) and four values for d 2 n{< 2 fm) with vls = 0 MeV. 
Then, the number of bases used to describe '^"Be in the 
present full calculation is (8 -I- 4) x 3 x 3 = 108. 

In this study, we would like to examine the effect of 
core structure on dineutron correlation in '*^Be. To this 
end, we compare the results obtained by two types of 
calculations. One is the full calculation mentioned above 
where the bases with da = 2, 3, 4 fm are superposed to 
describe the fluctuation of the a-a distance in the core 
structure, i.e., the shape fluctuation of the core. We de¬ 
note the full calculation as the “full-da” calculation. The 
other is the truncated calculation where only the bases 
with da = 2 fm are superposed, which gives the state 
with the smaller core deformation. We denote this cal¬ 
culation as the “fixed-da” calculation. The bases in the 
fixed-da calculation are recalculated because we use the 
interaction parameters modified from those of the full- 
da calculation to reproduce the two-neutron separation 
energy in each calculation. 

D. 2n overlap function 

For analysis to examine the dineutron correlation 
around the core, we propose a method using a 2n overlap 
function. The 2n overlap function / is defined as 

fir, re) = A{<^is,J{rA-i - VnJxtTn 

xd(rA - r„Jx4.T-„ |$ioBe(o-H)), (11) 

r = - T-ni, re = (r„i-f r„J/2, (12) 

where is the ®Be core wave function (Eq. dS])) spec¬ 
ified by the parameter da- We use the core with da = i 
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fm in the full-do, calculation and that with do = 2 fm 
in the fixed-do calculation, which gives the largest over¬ 
lap with the ground state obtained in each calculation, 
/(r, ro) is defined as a function of the relative and center 
of mass coordinates of two neutrons defined by the two- 
neutron coordinates, and , as shown in Eq. ©■ 
We omit the recoil of the core and locate its center of 
mass in the bra and ket states at the origin to measure 
the largest overlap between these states. The quantity 
f{r,rc) brings out information about the spatial distri¬ 
bution of spin-up and -down neutrons in <i>ioBe as a func¬ 
tion of r and rc- This quantity corresponds to a type 
of reduced width amplitude extended to the three-body 
case. Although the definition of the two-body density 
is nontrivial, we define /(r, rc) to be totally antisym¬ 
metrized, and f{r,rc) directly reflects the information 
about the spatial distribution of the two neutrons. The 
2n overlap function is useful for analysis of the dineutron 
correlation, as accomplished using the three-body model. 

In particular, we are interested in the component of 
the dineutron where two neutrons are coupled to a spin 
singlet with a relative s wave {I = 0). We focus on the 
dineutron in the S wave {L = 0) with respect to the core 
for simplicity. Here, we use the label “s” {I = 0) for 
the angular momentum for the relative coordinate r and 
“S” (L = 0) for the angular momentum for the center 
of mass coordinate tq. Thus, we project the angular 
momenta of the relative and center of mass motions of 
the two neutrons to the s and S waves, respectively, and 
consider the 2n overlap function, for the 

states projected to Z = L = 0 as a function of r = |r| and 
fG = kd- The details of the calculation of are 

given in the Appendix. 

Using the 2n overlap function calculate the 

probability 

Pf=rL=o = J r^drrldrc l/zlliokuc) |^ (13) 

which corresponds to the component of the spin-singlet 
and relative s-wave 2n pair moving in the S wave around 
the core. 

We also calculate the root-mean-square distance for r 
and rc of the two neutrons in the 5" = 0 and I = L = 0 
components: 

\/(-df=t=o 

= (^J ddrr^drGd|/fjl%(r,rG)r/^df=o) ' , 

(14) 

\J {'^‘g)i=l=o 

(15) 

which correspond to the dineutron size and the expansion 
from the core, respectively. 


Note that the asymmetry with respect to the exchange, 
r/2 -O' TG in reflects the mixing of single¬ 

particle orbits for two neutrons of different parity in the 
S' = 0 and I = L = 0 components. We decompose 
and HjfdLo fol® symmetric and antisymmetric compo¬ 
nents as follows: 


fitL=o{r, ro) 

= ro) + /zll=o “V, ro), ( 16 ) 

pS—0 P'S'—0(++) I pS—0( ) 

where //t(L=o Tz=z,=o defined as 


(17) 




= 9 (/ilL=o(u?'G) ± fitL=oi‘^rG,r/2)) , 

( 18 ) 


P, 


,S=0(±±) 


l=L=0 


= r^dr r^drc 


fS=0{±±) 

h=L=0 


{r, ra) 


(19) 


Here 'Pr/2-H-rG ^ Ih® ?'G exchange operator and 

the projection (1 ± 7^r/2<->T-G )/2 is equivalent to the dou¬ 
ble projection (1 ± Vr„-^)l‘2 x (1 ± Pr„ 2 )l‘^ H^® single¬ 

particle parities on /;di!'=o(^/2, ?'g), as shown in the 
Appendix. 7^r„i,2 are the space reflection operators of 

’’ni.s This means that (/ilk=o”^) 

and (-Pi=z,=o indicate the contributions of 

pure positive- (negative-) parity single-particle states of 
two neutrons in and We hereafter label 

/;=”l=o as the (-1—1-) or (-) component of the 2n over¬ 

lap function. The dineutron correlation in ^*^Be is seen 
in as the coherent mixing of the minor 

into the major \ as shown in Sec. IHI Cl 


III. RESULT 
A. Effective Hamiltonian 

In the present work, we use the Hamiltonian 

= T-rG + foce„t + ULS + foCoul, (20) 

where T and Tq are the total and center of mass ki¬ 
netic energies. In the present framework, since all the 
Gaussian widths in the present AMD wave functions are 
not equal, the center of mass motion cannot be removed 
exactly. We therefore treat the center of mass motion 
approximately by reducing the expectation value of Tq 
from the total Hamiltonian. Vcoui is the Coulomb force 
that is approximated by the summation of seven Gaus- 
sians. Ucent and Vls are the effective central and spin- 
orbit interactions. We use the Volkov No.2 force as 
14ent and the spin-orbit part of the GSRS force as 
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TABLE I: The used Majorana parameter m, the two neutron 
separation energy S 2 n, the root-mean-square radii rm,p,n of 
matter, protons and neutrons, and the expectation value of 
the squared neutron total spin (S^) in the full-do, (full) and 
fixed-do (fixed) calculations. The experimental values of the 
matter radii are referred from Ref. 1^. 



m 

S 2 n (MeV) 

rm (fm) 

rp (fm) 

c„ (fm 

0 (5^) 

full 

0.60 

6.71 

2.43 

2.22 

2.51 

0.39 

fixed 

0.64 

6.34 

2.33 

2.01 

2.38 

0.50 

Expt. 


8.48 

2.30 ±0.02 





In this work, we choose the strength of the spin- 
orbit force to be uls = 1600 MeV, as has been used in 
the previous works on the subject of ^’’Be [dllld^. The 
Bartlett and Heisenberg parameters in the central force 
are b = h = 0.125 which reproduce the deuteron binding 
energy and the n-n unbound feature. We use the Majo¬ 
rana parameter to = 0.60 as it was used in Refs. 
in the full-da calculation. In the fixed-da calculation, we 
choose TO = 0.64 to give almost the same two-neutron 
separation energy, 5'2n, that was obtained in the full-da 
calculation. 

B. fundamental properties 

We calculate the two-neutron separation energy, S 2 n, 
the matter, proton and neutron radii, rm,p,n, and the ex¬ 
pectation value of the squared neutron total spin, (5"^), 
shown in Table IH The two-neutron separation energy is 
calculated as the difference between the total binding en¬ 
ergy and the core binding energy, 

S2n = -((‘I’iOBe(0-|-)|I^I‘I’iOBe(0+)) 



where $8Be(o-i-) is the superposition of the ®Be core wave 
functions with da = 2, 3, 4 fm in the full-da calculation 
and that with da = 2 fm in the fixed-da calculation, and 
they are projected onto = 0^. The binding energy 
of ^°Be in each calculation is —60.42 MeV (full-da) and 
—49.63 MeV (hxed-da). The root-mean-square radii of 
matter, protons and neutrons are larger in the full-da 
calculation than those in the fixed-da calculation. This 
is natural because the core size becomes larger in the full- 
da calculation due to the fluctuation in the a-a distance. 

We show the neutron spin expectation value, in 

Table m In the present calculation, the ®Be core has zero 
proton- and neutron-spins and, therefore, (5"^) indicates 
the squared spin expectation value of the two valence 
neutrons. The finite value of (S”^) reflects the spin-triplet 
component; in other words, the degree of the LS disso¬ 
ciation of a dineutron. iS'^ is larger in the fixed-da cal¬ 
culation than that in the full-da calculation, indicating 
that the LS dissociation increases when the core struc¬ 
ture is fixed to be small. This point is discussed later in 


connection with the dineutron enhancement due to the 
core structure change. 

We have checked the relationship between the present 
AMD wave function and the DC wave function used in 
our previous work [d^. For the DC wave function, we 
assume a spin-singlet 2n pair around a core, and we 
have superposed the DC wave functions with the AMD 
wave functions. The main role of the DC wave function 
discussed in Ref. [i^ is to describe the dineutron-tail 
component in ^®Be. If we superpose the DC wave func¬ 
tions used in Ref. with the present AMD wave func¬ 
tions, the ground state energy is lowered by only 300 
keV at most and the other properties are largely un¬ 
changed. This means that the contribution of the DC 
wave function is minor, and that the structure as well as 
the dineutron tail can be well-described with the present 
(i 2 n-constrained AMD wave functions, at least for a ^°Be 
system with sufficiently bound valence neutrons. 


C. Effect of core structure change on dineutron 
correlation in ^°Be 

In this section, we discuss the dineutron correlation in 
the ground state of ^'^Be, focusing mainly on the effect of 
®Be core structure change. 

As an example of the analysis of dineutron correla¬ 
tion using the 2n overlap function, we plot the 2n over¬ 
lap functions, full-da cal¬ 

culation as functions of (r/2,rG) in Fig. [1] To show 
the asymmetry with respect to the tq = c/2 line, we 
show the r 12 -tg plot instead of the r-rc plot. There are 
two peaks in /jfr/JLo /;=”l=o ^ three peaks in 

/;=”l=o ■ Hereafter, we refer to the peak in the tq > r 12 

region of as the dineutron peak and the peak in 

the re < c/2 region of as the cigar peak. The 

two-peak structure in ^ comes from the two neu¬ 

trons occupying the orbits and the three-peak struc¬ 
ture in comes from those occupying the (sd)^ 

orbits. The (-) and (++) components of the 2n over¬ 
lap function, ^ and are antisymmetric 

and symmetric with respect to the re = c/2 line, re¬ 
spectively. In general, without the mixing of different 
parity contributions for the single-particle orbits, the ab¬ 
solute amplitudes of the dineutron and cigar peaks are 
exactly the same, and we do not describe the case as 
dineutron correlation. In other words, dineutron correla¬ 
tion is reflected in the asymmetry between the dineutron 
and cigar peaks with respect to the re = c/2 line in 
In the present case, the asymmetry in /;f!2*Lo 
comes from the mixing of the minor (sd)^ component 

(/j=”l=o) iiifo the major component as 

shown in Eq. (HSl). The relative phases at the dineutron 
peaks (c/2 ^ 1 fm) in and are coher¬ 

ent; on the other hand, those at the cigar peaks (cg ^ 1 
fm) in and are incoherent, leading to 
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FIG. 1: (Color online) The 2n overlap functions of (a)/;'L^'Lg, 
(b)/j‘^/jJ^Q and in the full-dc calculation as 

functions of (r/2,rG). For the guide, the colors of the lines 
corresponding to the region where the absolute amplitude is 
largest, intermediate and smallest are black, gray, and white, 
respectively. 


asymmetry between the dineutron and cigar peaks. The 
dineutron enhancement is seen in the obvious asymme¬ 
try between the dineutron and cigar components com¬ 
ing from significant mixing of the different-parity single¬ 
particle orbits for each valence neutron. 


We now investigate the dependence of the degree of 
dineutron enhancement on the core structure by com¬ 
paring the results obtained from the full-da and fixed-da 
calculations. In Table. im we show some properties of 
the 2n pair in two calculations. As mentioned above, the 
mixing of the different-parity single-particle orbits for the 
two neutrons reflects the dineutron correlation, in other 
words, the mixing ratio of into ^ re¬ 

flects the strength of dineutron correlation. The mixing 

ratio ^V^i=L=o ^ full-da calcu¬ 

lation and 3.0 % in the fixed-da calculation, meaning 
that the large deformation and large shape fluctuation 
in the ®Be core structure enhance the dineutron corre¬ 
lation. The stronger dineutron correlation in the full-da 
calculation is also reflected in the slightly smaller value 

of \/{^'^)f=L=o tvs'll that in the fixed-da calculation. 

The reason for which the core structure change en¬ 
hances the dineutron correlation is as follows; when the 
core deformation is fixed to be small in the fixed-da calcu¬ 
lation, two valence neutrons are distributed near the core 
to a larger extent and they feel the stronger spin-orbit po¬ 
tential at the surface. In the ^°Be case, the TS'-favored 
orbit of 0 p 3/2 is partially unoccupied and the valence neu¬ 
trons are favored energetically to occupy the 0 p 3/2 orbit. 
Under the smaller core deformation, the simplest shell- 
model component (two neutrons occupy only the lowest 
shell) is predominant because of the spin-orbit potential 
from the core and, therefore, dineutron correlation is not 
greatly enhanced because of the LS dissociation effect. 
On the other hand, when the core deformation becomes 
larger, the mean field generated by the core expands to 
the farther region so that two valence neutrons can be 
radially expanded from the core to form a spin-singlet 

compact dineutron. As a result, ^{rQ)fP^^Q is larger in 
the full-da calculation. Moreover, the spin-orbit poten¬ 
tial becomes weaker in the region far from the core, re¬ 
sulting in the suppression of the LS dissociation effect on 
the dineutron. The suppression of the LS dissociation in 
the full-da calculation is seen by the smaller spin-triplet 
component {{S^)/2) than that in the fixed-da calcula¬ 
tion, as already shown in Table IH 

However, it should be noted that the value 

itself is smaller in the full-da calculation than that in 
the fixed-da calculation. This is due to the fact that, if 
the core deformation becomes larger, the spin-singlet 2n 
components in the strong-coupling channels between the 
core and 2n (e.g. (Tcore = 2) 0 (A = 2) = 0 which is pro¬ 
jected out in the calculations of Pi^^-q and 
are mixed to a greater extent. We emphasize that the 
mixing of the (-1—h) component of the single-particle or¬ 
bits of 2n into the dominant (-) component increases, 

reflecting the enhancement of the dineutron correlation. 

To investigate the effect of the core structure on the 
2n properties in more detail, we compare the 2n over¬ 
lap functions in the full-da and fixed-da calculations. In 
Fig. m we plot the absolute values of the 2n overlap 
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TABLE II: The probabilities of the spin-singlet 2n pair those of the positive- or negative-parity components 

and the root-mean-square expectation values of the relative and center of mass distances of the 2n pair and 

in the full-da (full) and fixed-da (fixed) calculations. 



pS-0 

-n=L=o 

pS—0{ ) 

M=L=0 

p5—0(++) 

M=L=0 

(fm) 

VirDi^Lo (fm) 

full 

0.515 

0.491 

0.025 

3.70 

2.53 

fixed 

0.546 

0.530 

0.016 

3.75 

2.37 



0 1 2 3 4 5 6 

r/2 (fm) 


FIG. 2: (Color online) The absolute values of the projection of 
the 2n overlap functions of onto the r/2 axis in the full- 

da (full) and fixed-da (fixed) calculations. The red solid lines 
correspond to the full-da calculation and the blue dashed lines 
correspond to the fixed-da calculation. The thick lines are the 
dineutron component (the amplitudes in the r/2 < re region) 
and the thin lines are the cigar component (the amplitudes in 
the r/2 > ro region). 

function, |/;f!2*Lol’ projected onto the r/2 axis obtained 
in the full-da and fixed-da calculations. It can be seen 
that the dineutron peak (r/2 ^ 1 fm) is larger than the 
cigar peak (r/2 ^ 2 fm) in both calculations. However, 
the difference between the dineutron and cigar peaks is 
larger in the full-da calculation than that in the fixed-da 
calculation. This means that the dineutron correlation 
is enhanced due to the fluctuation in the a-a distance, 
as discussed in connection with the mixing of 
into It should be noticed that the absolute 

amplitudes of both the dineutron and cigar peaks them¬ 
selves become smaller in the full-da calculation because 
of the mixing of the (Lcore 0) 0 (L 0) = 0 compo¬ 
nents, as mentioned above. Introducing the fluctuation 
in the distance between 2q:s, the cigar peak decreases 
more than the dineutron peak and the difference between 
these peaks is certainly increased, indicating dineutron 
enhancement depending on the core structure change. 

D. Utility of the d 2 n-constraint calculation 

Finally, to show the utility of the present d 2 n- 
constraint calculation, we compare the dineutron-tail 
component obtained by the full-da calculation with that 



0 1 2 3 4 5 6 

rg (fm) 


FIG. 3: (Golor online) The absolute values of the projection 
of the 2n overlap functions of onto the ra axis in the 

full-da calculation and the calculation using the basis that 
gives the minimum energy. The red solid lines corresponds 
to the full-da (full) calculation and the green dashed lines are 
the calculation using only the minimum-energy basis (min.). 
The thick lines are the dineutron component (the amplitudes 
in the r/2 < re region) and the thin lines are the cigar com¬ 
ponent (the amplitudes in the r/2 > ro region). 


obtained by the single basis having the parameters 
{da,d 2 n,V 2 n) = (3,1.8,0.235), which is the minimum- 
energy AMD wave function with the assumption of a 2a 
core. The minimum-energy basis has an overlap of ~ 82 
% with the state obtained by the full bases calculation, 
and its energy is —55.80 MeV, which is ^ 5 MeV higher 
than the ground-state energy in the full-da calculation. 
This indicates that the superposition of the basis wave 
functions significantly improves the ^°Be wave function. 
We plot in Fig. [3] the absolute values of the 2n overlap 
function projected onto the tq axis for the states 

obtained in the full-da calculation and the minimum- 
energy basis calculation. In the calculation using the 
minimum-energy basis, the amplitudes of the dineutron 
{rc ~ 2 fm) and cigar {tq ~ 1 fm) peaks do not differ 
greatly, because the 0 p 3/2 orbit (the lowest shell-model 
limit) is dominant for 2n and the dineutron correlation 
mostly disappears in this basis. 

A remarkable difference between these calculations is 
seen in the dineutron-tail component {rc ^ 4 fm). The 
minimum-energy basis does not have the dineutron tail 
expanded into the farther region because in a single ba¬ 
sis, this tail has a rapidly dumping Gaussian form. On 
the other hand, in the full-da calculation where many 


















bases with various d 2 n and V 2 n values are superposed, 
the dineutron tail is improved remarkably. This indicates 
that the present framework well-describes the dineutron 
tail. It is also expected that this method can be useful 
for extremely loosely bound nuclei with neutron-halo or 
-skin structures. Further application to loosely bound 
nuclei is future work. 


IV. SUMMARY 

In this study, we extended the AMD framework to de¬ 
scribe dineutron correlation around a core with various 
structures. We first prepared the core wave functions 
using the AMD method, which was useful for describ¬ 
ing various structure changes such as deformation and 
clustering. We described the detailed two-neutron mo¬ 
tion around the cores considering the relative distance 
between the core and the two neutrons as the degree 
of freedom, and we superposed the basis wave functions 
with various core-2n distances to describe the competi¬ 
tion between the dineutron formation and the LS dis¬ 
sociation at the nuclear surface as well as the dineutron 
expansion far from the core. We additionally changed 
the Gaussian widths of the two neutrons to take into 
account the dineutron dispersion effect. Moreover, we 
have constructed a 2n overlap function as the analyzing 
method for investigating dineutron correlation around a 
core in detail. This method enabled us to visualize the 
spatial correlation between two neutrons and made the 
discussion clear. 

As a first step, we applied the framework to ^°Be and 
focused on the effect of the ®Be core structure change 
on the dineutron correlation. In the present work, we 
assumed a 2a cluster structure for the ®Be core and con¬ 
sidered the fluctuation in the a-a distance as the core 
structure change. Dineutron correlation was seen in the 
degree of mixing of different-parity single-particle orbits, 
i.e., the mixing of the minor (sd)^ component into the 
major component for two neutrons in the ^*^Be case. 
When the core was less deformed, the two neutrons were 
distributed near the core and greatly dissociated to the 
independent ( 0 ^ 3 / 2 )^ component due to the spin-orbit 
potential from the core; as a result, dineutron correla¬ 
tion was suppressed. On the other hand, when the core 
was well-deformed by taking into account the fluctua¬ 
tion in the distance between 2QfS, two neutrons could be 
expanded radially and they were favored to form a dineu¬ 
tron. Moreover, at the region far from the core, the spin- 


orbit potential becomes weaker and the LS dissociation 
effect was suppressed. As a result, the dineutron correla¬ 
tion was enhanced at the surface and further regions due 
to the core structure change. 

In this work, we have applied the extended methods 
to ^°Be and showed that the present framework is useful 
to describe the dineutron component around the well- 
deformed core. Additionally, we have concluded that the 
core structure significantly affects the dineutron correla¬ 
tion. Certainly, ^°Be is not a very loosely bound nucleus 
and the LS dissociation effect on the dineutron plays 
an important role. However, we could see the obvious 
asymmetry between the dineutron and cigar components 
(Fig. [1]), and would like to stress that dineutron correla¬ 
tion can be more or less apparent in most neutron-rich 
nuclei, even in those without an outstanding neutron- 
halo or -skin structure. In the future, we will apply the 
present framework to various nuclei and clarify the uni¬ 
versal properties of dineutron correlation, e.g., the effect 
of the core excitation and clustering in addition to that of 
the core deformation on dineutron correlation, through 
the systematic investigation. 


Appendix: Overlap functions / and 

Here we show the details of the overlap functions / and 
j(±±). We define the general 2N overlap function / as 
below: 

f{r,rG) = M($core<5(T'A-l “ TnJxA-ITA-I 

xS{rA - rN2)XATA (A.l) 

r = rjVa - rN^, re = {tni + r'Ar2)/2, (A.2) 

where $totai is the total wave function with A nucleons 
whose center of mass of the core component is located 
at the origin and <I>core is the core wave function with 
A — 2 nucleons whose center of mass is located at the 
origin. XA-i,A and ta-i,a are the spin and isospin wave 
functions of the {A — I)th or Ath nucleon and we can 
choose t or 4 , for the spin components and p or n for 
the isospin components for each purpose. In the present 
work, we want to investigate the dineutron motion that 
is a spin-singlet pair of two neutrons so that XA-i =t 
and XA =4-, and ta-i,a = n. In this appendix, we show 
only this case but generalization is simple. 

In the case of a spin-singlet 2n pair, Eq. (lA.lIl can be 
rewritten as 




(A.3) 


x(xtlx*)(X4.IXj) X 
/4 

= + ’'V4) - - i^jyf - {vi - Uj)r ■ rG/2] 

i,j^n ' ' 

- ’^jYjy)Yimir)Y;^{Yr) 

I m 

+ i'jYj)rG)Yim{rG)YCn{YG) 

I m 

x(xtlx*)(X4lXj) X det^P’^), (A.4) 


where j; are the spherical Bessel functions and Yim 
are the spherical harmonics. Yj-^g the polar an¬ 
gles of the vectors UiYi — VjYj or ViYi + I'jYj, re¬ 
spectively. is an (A — 2) x (A — 2) norm matrix 

composed of (v^core^l (k = 1,..., A - 2) in (<i)core| and 
|</3totai,K) (k = 1,..., A) in |$totai) except for the Ah and 
jth single-particle wave functions, /(r, re) is totally an¬ 


tisymmetrized by A and hence yvh per¬ 

form the angular integrals J (PrYQf^{r) f (Pi'GYooif'G) to 
project the relative and center of mass motions of 2n onto 
I = L = 0, and we additionally neglect the term propor¬ 
tional to r • re in the exponential term in Eq. (IA.4I) . 
resulting in /ilA°=o(A ^ 0 ) as 


J 


/ili=o(A tg) = + y/y - '^iY? - VjY^] 

i,j^n ^ 

x47rjo {-i{iyiYi - VjYj)r) x jo {-2i{viYi + VjYj)rG) 

x(xtlx*)(XilXi) X detSP’^). (A.5) 


This 2n overlap function only depends on the abso¬ 
lute values of r = |r| and tg = vg- We calcu¬ 
late the probability of the spin-singlet 2n pair, 
the root mean square distance between two neutrons 

1 =^= 0 ’ between the two neutrons and the 

core (^g)£l*=o simply using the numerical integrals of 
r and rg (Eqs. (IT51) . (fTTl) and (fTCll l. 

We calculate the (++) and (-) components 

/i=”L=o^^ from the 2n overlap function as 


Eq. (fT^ . Here, we show that defined in 

Eq. (1181) . is certainly equivalent to the components where 
both neutrons are projected to positive- or negative- 
parity single-particle orbits if the relative and center of 
mass motions of 2n are projected as / = L = 0. We 
begin by projecting both two neutrons in the 2n over¬ 
lap function f{r,rG) to the positive- or negative-parity 
single-particle orbits. Noting the definition of the coor¬ 
dinates r and tg in Eq. (IA.2I) . we have 


(1 ± J /2 X (1 ± J /2 X/(r, re) 

= ^ [(1 + 'Pr-^.'Pr^y ± {rr„, + f {v, Tg) 

= ^ lify, rG) + fi-r, -tg)) ± (/(2rG, r/2) -h /(-2rG, -r/2))\. (A.6) 


If we project f to I = L = 0 with the operator Vi=l=o, 


becomes just a function of the absolute values of 
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r and tq {'Pi=L=of{r,rG) = As a result 

of the projection to / = L = 0, the former two terms and 
the latter two terms on the last line of Eq. (jA.6D give 
the same contributions, respectively. We can therefore 
rewrite the Eq. (IA.6I) under the projection to Z = L = 0 
as 

ri^L=o [{l±rr„J/2 X (1 ±7’.„J/2 X f{r,rG)] 

= ^ iff=L=oir,rG) ± /ill=o(2rG,r/2)) 

= (A.7) 

This means that defined in Eq. (fl^ are cer¬ 

tainly the components in the 2n overlap function where 


both two neutrons are projected to the positive- or 
negative-parity single-particle orbits if the relative and 
center of mass motions of the two neutrons are projected 
as I = L = 0. 
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